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Particle production of identified charged hadrons, tt*, K"^, p, and p in Au+Au collisions at 
-y/ijvjv ~ 200 GeV has been studied as a function of transverse momentum and collision centrality 
at y = and y ~ 1 by the BRAHMS experiment at RHIC. Significant collective transverse flow 
at kinetic freeze-out has been observed in the collisions. The magnitude of the flow rises with the 
collision centrality. Proton and kaon yields relative to the pion production increase strongly as the 
transverse momentum increases and also increase with centrality. Particle yields per participant 
nucleon show a weak dependence on the centrality for all particle species. Hadron production 
remains relatively constant within one unit around midrapidity in Au+Au collisions at ^/sNN = 
200 GeV. 

PACS numbers: 25.75.Dw 
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I. INTRODUCTION 

The primary goal of the relativistic heavy-ion collider 
(RHIC) is to create and study matter at extremely high 
energy density. It is hypothesized that at the energy den- 
sities reached in central Au+Au reactions at RHIC, the 
matter created is composed of de-confined colored ob- 
jects 0, 01 13 ■ A summary of the results and opinions 
of the four experimental collaborations on the status of 
achieving this goal can be found in their "White Papers" 
0, IE 111, 13] • We expect that the signals of any de-confined 
phase should become stronger as the overlapped region 
increases in Au+Au collisions. Testing this hypothesis 
requires studying particle production as a function of 
centrality. The particle distributions in transverse mo- 
mentum and rapidity may provide a key to understand- 
ing any non-hadronic effects that might appear in central 
nucleus-nucleus collisions. Pions, kaons, protons and an- 
tiprotons are the most abundantly produced particles in 
the high-energy heavy-ion collisions, and they carry in- 



formation about the bulk properties of the nuclear matter 
created from the collisions. 

Pions, being the hghtest of the produced hadrons are 
thus the most copiously produced, and their numbers can 
be related to the entropy density of the emitting source. 
Kaons carry a significant fraction of the total strangeness 
produced Ig]- Protons and antiprotons provide an ex- 
perimental tool for measuring baryon production and al- 
low us to explore baryon transport from beam rapidity 
toward midrapidity . The global thermodynamic 

properties and collective motion of the system at the ki- 
netic freeze-out point can be deduced, albeit, in a model 
dependent way, from transverse momentum spectra as a 
function of rapidity and centrality 

In this paper, we present transverse momentum spec- 
tra, yields, and ratios for identified charged hadrons (tt^, 
K^, p, p) obtained with the BRAHMS Mid- Rapidity 
Spectrometer. We have measured these spectra for y = 
and y ^ 1 as a function of collision centrality. At midra- 
pidity, our observations are in agreement with the result 
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of the PHENIX experiment within systematic uncer- 
tainties. The data presented here are available at [T^. 

II. EXPERIMENTAL DETAILS 

The BRAHMS experiment consists of two movable 
magnetic spectrometer arms, the Mid-Rapidity Spec- 
trometer (MRS) and the Forward Spectrometer, and 
global detectors for event characterization. 

In order to characterize the centrality of collisions, a 
multiplicity array (MA) consisting of a coaxial arrange- 
ment of Si strip detectors and scintillator tiles surround- 
ing the intersection region is employed. The Si strip de- 
tectors and scintillator tiles give independent measure- 
ment of charged-particle multiplicities allowing the two 
measurements to be averaged in the final determina- 
tion. The pseudo-rapidity coverage of the MA is approx- 
imately -2.2<?7<2.2 The centrality selection 
is obtained by developing minimum-bias trigger events, 
which are defined using two Zero Degree Calorimeters 
(ZDC), requiring energy deposit equivalent to at least 
one neutron in each of the two detectors and also re- 
quiring a signal in the MA to reject Coulomb dissocia- 
tion events |l4l.ll6|. Figure shows the charged-particle 
multiplicity distribution for minimum-bias events in the 
range of— 2.2 < 77 < 2.2. The lines on the plot in- 
dicate four centrality windows in the analysis, 0—10%, 
10-20%, 20-40%, and 40-60%, where 0% corresponds 
to the most central events. 




Normalized Multiplicity 



FIG. 1: MA array multiplicity distribution. Lines show the 
limits for indicated centralities. 

The number of participating nucleons (Npart) in Ta- 
ble m are estimated using the Glauber Monte-Carlo HI- 
JING calculation T?]. More peripheral collision events 
were not included in this paper because of limited statis- 
tics. 

The uncertainty in determining centrality from the 
multiplicity distribution was estimated to be ±1.7% for 
the most central bin and ±9.4% for the most periph- 
eral bin. The fraction of the inclusive yield lost by the 
minimum-bias trigger is estimated to be about 4% and 



Centrality 


(Npart) 


- 10% 


328 ± 6 


10 - 20% 


239 ± 10 


20 - 40% 


140 ± 11 


40 - 60% 


62 ± 10 



TABLE I: The number of participant nucleons Npart esti- 
mated from the HIJING model jl^l- 



is corrected for. The location of the collision vertex was 
determined to an accuracy of 0.7 cm using Beam-Beam 
Counters (BBC) O. The BBCs are located 2.2 m on 
either side of the nominal interaction point (IP) and also 
provide a start time (TO) for time-of- flight measurement. 

The MRS is a single-dipole-magnet spectrometer 
which, by rotating about the nominal collision point, pro- 
vides the angular coverage of 30° < 9 < 95°, where 6 is 
the polar angle with respect to the beam axis. The MRS 
contains two time projection chambers (TPCs), TPMl 
and TPM2, which determine the three dimensional tra- 
jectories of the charged particles through the spectrom- 
eter. Between the two TPCs there is a dipole mag- 
net (D5) for momentum determination. This assembly is 
followed by a highly segmented scintillator time-of-flight 
wall (TOFW). 

The BRAHMS TPCs are located at a distance 
0.95 m (TPMl) and 2.85 m (TPM2) from the interaction 
point. Each TPC is a rectangular box filled with 90% Ar 
and 10% CO2. The ionization produced by charged par- 
ticles is collected on an anode grid. This grid is divided 
along the particle path into 12 rows (TPMl) and 20 rows 
(TPM2). Each row has 96 pads (TPMl) and 144 pads 
(TPM2) transverse to the direction of a normal-incident 
particle. The mapping of row, pad, and drift time leads 
to three-dimensional space points. The averaged resolu- 
tions measured from track residuals are 310—387 /im for 
X (pad) and 427— 490_um for Y (time). Details can be 
found in Ref. [Hill El 13. 

Track reconstruction starts by finding straight-line 
track segments in the TPCs. The reconstructed straight 
tracks are joined inside the analyzing magnet by taking 
an effective edge approximation, and the momentum as- 
sociated with the tracks are calculated from the vertical 
magnet field, the length in the magnetic field region, the 
polar angle of the tracks with respect to the matching 
plane, and the averaged vertical slope of the tracks. The 
matching plane is defined as the vertical plane that con- 
tains the perpendicular bisector of the line joining the ef- 
fective edge entry and exit points of the tracks. Once the 
momentum is known the reconstructed tracks are pro- 
jected toward the beam axis and checked for consistency 
with the collision vertex determined using the BBCs. For 
this analysis we only use tracks that project to within 
±12.5 cm of the nominal vertex in the horizontal plane. 

The momentum resolution of the spectrometer can be 
inferred from the width of our mass-squared distribu- 
tions since — p'^[{ct/l)^ — 1], where c is the velocity 
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momentum [GeV/c] 

FIG. 2: Momentum resolution 5p/p calculated for the MRS 
90° setting as a function of momentum with a magnetic field 
of 6 kG. The curves are obtained from parametrization of the 
mass resolution o"^2- At low momentum multiple scattering 
dominates the resolution while above 1 GeV the angular res- 
olution of the spectrometers is the most important effect. 



of light, I the track path length, and t is the time-of- 
flight whose resolution is measured independently. The 
extracted resolutions were fit to the form (Sp/p)^ = 
(cresP)'^ + (cmuiti/ P)'^ , wherc Cres is the Contribution from 
the intrinsic angular resolution of the tracking detectors 
and Cmuiti is the resolution from multiple scattering. The 
best fit is given by Cres ~ 0.014 c/GeV and Cmuiu ~ 0.01 
with D5 at 6 kG. Figure |21 shows resolution curves based 
on this fit for pions, kaons and protons when the D5 mag- 
net is set to 6 kG. For the data presented in this paper 
the momentum resolution lies between 2 and 8% depend- 
ing on the momentum of the particle and the magnetic 
field in D5. 

Particle Identification (PID) is based on time of flight 
data from the TOFW {at ~ 80 ps) with the start time 
taken from the BBC {at » 35 ps). The TOFW is located 
at a distance of 4.3 m from the nominal interaction point. 
Final mass identification PID is based on cuts in the 
vs. p space, as shown in Fig. |2| The cut boundaries are 
set at ±3tT from the mean values. TOFW provides 
tt/K separation to the momenta of 1.85 GeV/c and K/p 
separation to the momenta of 2.85 GeV/c. For parti- 
cle identification above these momenta, the yields that 
are within the ±3(t overlapped regions are corrected for 
misidentified particles up to the momentum of 2 GcV/c 
for tt/K, and the momentum of 3 GeV/c for K/p. The 
kaon contaminations within the pion sample at 2 GeV/c 
and within the proton sample at 3 GeV/c are less than 
1%. The dotted fines in Fig. 13 show the upper momen- 
tum limit used in the analysis on each particle type. 

The data presented here were collected with the MRS 
at 90°, for y ~ (-0.1 < y < 0.1) and at 45° - 35° for 
y - 1 (0.7 < y < 1.1) in RHIC Run II (2001-2002). 

From the number of identified particles, invariant dif- 
ferential yields are obtained from several spectrometer 
settings for different collision centralities. The invari- 




Momentum [GeV/c] 

FIG. 3: Mass-squared distribution in TOFW as a function of 
reconstructed particle momentum obtained from tracking in 
TPCs through the D5 magnet at y = 0. The curves in the 
figure are PID boundaries set at ±3cr^2 from the mean of 
the particle. The dotted lines indicate the momentum cutoff 
for particle separation. 



ant yields are corrected for geometrical acceptance and 
the efficiencies for detecting particles in the spectrome- 
ter. The inefficiencies arise from two effects, inefficiencies 
due to single track loses and those due to multiplicity de- 
pendent effects. 

Geometrical acceptance factors are obtained from 
the GEANT ^ simulation package BRAG (BRAHMS 
Analysis Geant), which is based upon the geometry and 
tracking capabilities of the the BRAHMS experimental 
setup. The acceptance correction is calculated for each 
MRS setting and five different vertex windows covering 
the MRS track vertex range used in the analysis. 

The single track efficiency as a function of momen- 
tum in the spectrometer is determined by a Monte Carlo 
simulation. Events with one particle are first processed 
through BRAG with multiple scattering, decays, and 
hadronic interactions processes included. In order to 
evaluate these effects, the simulated events are processed 
through the same digitization, reconstruction and parti- 
cle selection algorithm that is applied to the real data. 

The upper panel of Fig. ^ shows correction factors ap- 
plied to the pion, kaon and proton spectra at y=0 to 
account for multiple scattering and (tt and K) decay in 
flight. The low momentum p spectra are also corrected 
for hadronic absorption in the beam pipe and detector 
materials. This effect amounts to ~ 2 — 3% of the to- 
tal yield. No corrections are applied for secondary pro- 
tons, arising mainly from the interaction of pions with the 
beam tube, since the contributions is found to be negli- 
gible in MC simulations using HIJING as input, when 
the tracks were required to point back to the IP. The 
difference of correction factors at different spectrometer 
settings is below 1%. 

The multiplicity dependent track reconstruction effi- 
ciency has been studied by embedding simulated tracks 
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inefficiencies. 

Protons and antiprotons from weak decays lead to 
a contamination of the primary hadron spectra. The 
proton and antiproton spectra are corrected to remove 
the feed down contributions from A and A weak de- 
cays. At midrapidity the ratio A^(A) = 0.89N{p), 
N{A) = 0.95N{p) has been reported in Au+Au coUi- 
sions at ys]^ 130 GeV HI El. We have studied 
the magnitude of the corrections using various model as- 
sumptions as input to the BRAG simulations. Assum- 
ing primary A/p ratios at y^SNN = 200 GeV similar to 
those measured at the lower energy 25] and a constant 
behavior with rapidity, we take these ratios and mea- 
sured spectra shapes as input to the BRAG code for feed 
down correction from A decays. The simulated tracks 
are generated for the full phase space, digitized, and go 
through the real data analysis algorithm, as is done to 
determine the other correction factors. The lower panel 
of Fig. 0] shows the ratio (e/eed) of secondary p and p to 
all measured p and |5 as a function of momentum. The 
fractional factors range from 25—40%, and the largest 
value is ~40% around a momentum p ~ 0.5 GeV/c. We 
multiply the proton and antiproton spectra by I— e feed 
for all centrality windows and rapidities as a function of 
momentum. The data are corrected on a track-by-track 
basis for efficiency and feed down contributions. 



FIG. 4: Correction factors applied to the invariant yields, 
as discussed in the text. The top panel shows the momen- 
tum dependence of the multiple scattering, decay in flight 
and hadronic absorption efi'ects. The bottom panel shows the 
fractional factors of p and p from A and A in all measured p 
and p, efeed, as a function of momentum. 



into real events at the raw data level j23|. The com- 
bined events are reanalyzed to determine if the embedded 
tracks are still reconstructed by the tracking program. 
Each track is associated with a number of pads in TPMl 
and TPM2. The resulting tracking efficiency is parame- 
terized as a function of the number of track-related pads 
found in the two TPCs with signals above threshold for 
pions, kaons, and protons in various spectrometer angle 
settings. The mean number of track-related pad hits in 
the data sample varies from ^350 to 60 as the central- 
ity varies from 0% to 60%. For the most central events, 
track reconstruction efffciency is ^85—95% depending on 
spectrometer angle setting. 

The efficiency for individual TOFW slats is investi- 
gated by projecting TPC tracks to the slat and compar- 
ing this to the distribution of TOFW hits. The possibility 
of having multiple hits on a single slat is also corrected 
for. The overall efficiency for particle identification using 
TOFW is estimated to be 90%±2% without significant 
dependences on spectrometer settings or collision central- 
ities. The inefficiency includes uncertainties in matching 
tracks with hits in TOFW as well as the intrinsic detector 



III. EXPERIMENTAL RESULTS 

Figure El shows the invariant spectra for charged 
hadrons, tt^, K^, p, and p at different collision central- 
ities. The overlaid lines indicate fits to the data from 
?/ = in the range shown. The pion spectra are fitted 
with a power-law function, j4(1 + Pt/po)~^ ■ For kaons 
the spectra are best fit by an an exponential in tut , 
^I'g-mT/Tj^ where niT — yPj^^rm^ where toq is the 
mass of the particle. The proton and antiproton rriT- 
spectra tend to deviate from a single exponential shape, 
so a sum of two exponential functions is used to in the 
fit. The point-to-point systematic uncertainties on the 
spectra and quality of fit are estimated by using other fit 
functions and varying the fit ranges. The errors shown 
on the data points are statistical only. The overall sys- 
tematic errors are estimated to be 10—15%. The main 
sources for the overall systematic errors are from the un- 
certainties in the normalizations used to calculate the 
invariant yields. Others are from uncertainties in es- 
timating background contribution, track reconstruction 
efficiencies, acceptance of spectrometer and particle iden- 
tification losses. The yields and mean transverse momen- 
tum values are extracted from the fit functions. Tables Hll 
gives the fit ranges and the estimated percentage of the 
total yield included in the fit ranges. 

Figure El shows the mean transverse momenta, (pr), 
as a function of Npart- We find that {pr) increases with 
particle mass and centrality. This is suggestive of a hy- 
drodynamic system where the initial pressure increases 
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FIG. 5: The invariant yield spectra for the identified par- 
ticles for the 0-10%, 10-20%, 20-40%, and 40-60% cen- 
tral collisions in Au-fAu collisions at y^ijvjv=200 GeV. Cir- 
cle(triangle) symbols are the results at y = (j/ ~ 1). Curves 
overlaid to the data points are fits to the y — data, as 
discussed in the text. For clarity, the data points are scaled 
down by a factor of 10 successively from the top (0—10%) to 
bottom(40— 60%) in decreasing order of centrality. The error 
bars are statistical only. 





y = 


y~ 1 


^± 


0.2 <PT< 1.9 (76.5%) 


0.3 <PT < 1.3 (72.1%) 




0.4 <pT < 1.9 (48.0%) 


0.3 <pT < 1.3 (40.9%) 


P,P 


0.3 <pT < 3.0 (72.1%) 


0.6 <pT < 2.1 (64.9%) 



TABLE II: Fit ranges for curves shown in Fig.|3 The yields 
were calculated from the data within the fit ranges, and the 
estimated percentage is the ratio of measured yields within 
the fit ranges to extrapolated yields for the full momentum 
range. 



with the number of participants. 

We also tried a hydrodynamic model fit to the exper- 
imental data with two free parameters: collective trans- 
verse flow velocity f3 and the global thermal freeze-out 
temperature Tf. We have utilized a version of a hydro- 
dynamically inspired "blast-wave" model initially devel- 
oped to describe lower energy data 0. Assuming ki- 
netic freeze-out of matter at constant Tf with a collective 
transverse flow characterized by a velocity /3s , the invari- 
ant rriT distribution can be described as follows: 
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FIG. 6: Mean pr extracted from the fit functions to the 
spectra as a function of centrality (Npart) for 7r+, , p (left) 
and 7r~, K~ , p (right). Open symbols represent at j/ = 0, and 
closed symbols are for y ^ 1. The error bars are statistical 
only. 



dN 



niT drriT 



r dr rriTloi 



PTsinhp m^coshp 



Tf ' Tf " 

(1) 

where Tf is the freeze-out temperature, /q, K\ are 
modified Bessel functions and p = tanh^^/Jx is the 
transverse rapidity. The transverse velocity profile /3t 
is parameterized by the surface velocity /3s : Pxir) = 
Ps{r/Rmax)°'- Results of a simultaneous fit for 0—10% 
centrality of Eq. ^ to the spectra for y = is shown 
in Fig. 13 The parameters a, Tf, and /3s were allowed 
to vary, as well as the normalization constants for each 
particle type. The px coverage of spectra at y ~ 1 is 
not sufficient for a reliable hydrodynamic model fit. The 
source parameter used was Rmax — 13 fm 26] . The inte- 
gral in Eq.^is relatively insensitive to changes in Rmax, 
changing by less than 5% when Rmax is changed from 5 
to 20 fm. 

For the most central events (0—10%), the fit yields val- 
ues of Tf = 109.6±0.9 MeV, /3s = 0.78±0.003, and a = 
0.40±0.05. The average flow velocity (/3t) is then esti- 
mated to be 0.65 by taking an average over the transverse 
geometry When fitting spectra from the other cen- 
trality windows we fixed the value of a to be 0.4. The 
fits in the four centrality bins give a x^/DOF between 
1.34—1.49. The systematic uncertainties in the fit pa- 
rameters are estimated to be less than 5%. Figure |H1 
shows the centrality dependence of the temperature and 
surface velocity. Tf decreases with centrality while /3s 
increases. Since the surface velocity keeps increasing un- 
til the system decouples these results suggest that central 
collisions decouple later. The increased energy associated 
with the surface velocity requires a lower final tempera- 
ture by energy conservation [28l| . 

The rapidity densities dN/dy are determined for each 
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FIG. 7: The parametrization and the rriT spectra for 0-10% 
centrality. The errors in the fit parameters are statistical only. 



y ■ 



- 10% 10 - 20% 20 - 40% 40 - 60% 



TT 

K+ 
K- 

P 

P 



283.3 ± 2.4 
277.9 ± 2.4 
45.0 ± 0.67 
40.9 ± 0.63 

18.6 ± 0.21 

13.7 ± 0.18 



196.2 ± 1.9 

195.3 ± 2.0 
29.8 ± 0.54 
28.1 ± 0.52 
13.4 ± 0.18 

9.7 ± 0.2 



119.7 ± 1.4 

118.8 ± 1.4 
17.9 ± 0.38 
16.3 ± 0.4 
7.6 ± 0.13 
5.6 ± 0.1 



47.5 ± 0.96 
46.3 ± 0.89 
6.3 ± 0.26 
5.9 ± 0.25 
2.81 ± 0.08 
2.22 ± 0.08 



- 10% 10 - 20% 20 - 40% 40 - 60% 



y • 



TT 

K+ 
K~ 

P 

P 



266.8 ± 2.4 
276.2 ± 2.7 
40.6 ± 0.32 
38.5 ± 0.33 

18.0 ± 0.17 

12.1 ± 0.14 



193.1 ± 2.2 
198.1 ± 2.5 
29.3 ± 0.27 

26.6 ± 0.29 

12.7 ± 0.15 
8.9 ± 0.13 



117.6 ± 1.3 
120.6 ± 1.5 
17.2 ± 0.17 
16.1 ± 0.19 

7.4 ± 0.1 

5.5 ± 0.09 



50.8 ± 1.01 
49.8 ± 1.0 
6.7 ± 0.12 
6.1 ± 0.13 
2.97 ± 0.08 
2.29 ± 0.08 



TABLE III: The yield dN/ dy from integration of extrapolated 
function in each centrality bin at y=0 and y ~ 1. The fit range 
was shown in Tables Ell The errors are statistical only. 



particle by integrating the measured px spectrum over 
Pt using the previously discussed functional form to 
extrapolate outside of the region of the measurement. 
Table IIIII shows the results at y = and y ~ 1, 
respectively. The dN/dy values per participant pair 
are shown as a function of Npart in Fig. |S| The sys- 



> 
O 



0.7 



FIG. 8: Kinetic freeze-out temperatures Tf and transverse 
flow velocities /3s resulting from a simultaneous fit of to tt^, 
, p and p spectra to Eq. Q as a function of collision cen- 
trality. The curves correspond to increasing centrality as /3s 
increases. The points represent the best fit values of Tj and 13s 
while the contours indicate Icr and 3(7 levels. The systematic 
errors are < 5% and are not included in the figure. 
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FIG. 9: Rapidity density (dN/dy) per participant pair 
[Npart /2) as a function of Npart for , , p, net pro- 
ton (left) and 7r~, K~ , p (right). Open symbols represent 
at J/ = 0, and closed symbols are for y ~ 1. The error bars 
represent the statistical errors only. 



tematic errors on the yields and (px) values from the 
extrapolation to the low momentum region are esti- 
mated as 5—10%. The only centrality dependence evi- 
dent is a small increase in the rapidity densities for the 
K and p channels in going to more central collisions. 
The net proton rapidity densities show a increase from 
40-60% to 20-40% and saturate after that. The proton 
excess,(A^p - Np)/{Np + Np), is 0.15±0.01-0.12±0.02 at 
y = and 0.19±0.01-0.13±0.02 at y ~ 1 from periph- 
eral to central collisions. Our p/p ratio from pp colli- 
sions showed a proton excess of 12% at midrapidity 2^. 
This baryon a.symmetry has been modeled at lower en- 
ergy systems ISOl where it has been found to be signifi- 
cantly greater [SlL Is^ . 
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FIG. 10; K'^ /tv'^ (left) and K~ /tv~ (right) ratios as a function 
of pt at 1/ = (top) and t/ ~ 1 (bottom). Closed and open 
symbols represent central (0—10%) and peripheral (40—60%) 
collisions, respectively. Only statistical error bars are shown. 



FIG. 11: p/n'^ (left) and p/tt~ (right) ratios as a function 
of pt at J/ = (top) and J/ ~ 1 (bottom). Closed and open 
symbols represent central (0—10%) and peripheral (40—60%) 
collisions, respectively. Only statistical error bars are shown. 



Figure [Till compares the K/tt ratios as a function of 
Pt for central and peripheral coUisions. The ratios for 
both charges increase with px at both rapidities. There 
is no significant centrahty dependence below 1 GeV/c, 
however, the increase is faster in central than periph- 
eral collisions at higher momentum. The weak centrality 
dependence of the K/tt production at RHIC [s^l differs 
from the measurements at lower energy from AGS [s^ l 
and SPS experiments j where the K~ /■k~ ratio is en- 
hanced 1.5 — 2) from central collisions compared to 
peripheral collisions. This might be attributed to the 
energy dep endent longitudinal geometry of the colliding 
nuclei |3(t1 IstII . but further experimental and theoretical 
work is needed to understand the observed energy depen- 
dence. 

Figure ITTI shows the p/7r+ and p/tt^ ratios as a func- 
tion of Pt obtained at y = and y ^ 1 for central and 
peripheral collisions. For both centralities the ratios rise 
fast at low Pt- Around pt ^ 2.0 GeV/c the ratios in- 
crease from peripheral to central collisions by 20% for 
p/iT~^ and by 50% for p/n~ ratios. Parton recombina- 
tion and quark coalescence models describe qualitatively 
the observed baryon to meson ratios for central colli- 
sions up through the intermediate pt region extending 
to ~4-5 GeV/c 



To clarify the centrality dependence of particle ratios 
at higher pT, and to see if the ratios of high pT particles 
are sensitive to the size of the interaction volume, we 
present the K/tt and p/n ratios versus Npart- Figure IT^ 
shows that the K/tt and p/ir ratios for 1.3 < pt < 2.0 
GeV/c increase with Npart- This increase is similar for 
protons and kaons with little difference in slope between 
the particle and antiparticle ratios. The behavior sug- 
gests that rescattering and/or hydrodynamic effects are 
stronger for larger collision volumes. 



IV. SUMMARY 

Particle production of identified charged hadrons, tt^, 
K^, p, and p in Au-j-Au collisions at y/s^N — 200 GeV 
has been studied as a function of transverse momentum 
and collision centrality at y = and y ~ 1 by the 
BRAHMS collaboration at RHIC. Significant collective 
transverse flow at kinetic freeze-out is estimated for both 
central and mid-central events. The magnitude of the ra- 
dial expansion increases with the collision centrality in- 
dicating more hydro-like collectivity in the transverse di- 
rection for the central collisions. The p, p and yields 
relative to the pion production at RHIC show a strong 
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FIG. 12: The ratio of K/n (left) andp/n (right) as a function 
of Npart for 1.3 < pt < 2.0 GeV/c at mid-rapidity. Open and 
closed symbols represent the ratio of positive particles and 
negative ones, respectively. 



transverse momentum dependence. Contrary to lower 
energy results no significant centrality dependent K/n 



ratios below 1 GeV/c are observed at the RHIC energy. 
The p/tt^ and p/tt^ ratios increase with px and central- 
ity reaching a value of 0.6 and 0.5 at ~ 2 GeV/c. The 
particle yield scaled by Npart is nearly constant, and only 
weakly increasing with centrality for all particles. No 
significant changes for the bulk properties in hadron pro- 
duction are observed within one unit around midrapidity 
in Au+Au collisions at y^SNN = 200 GeV. 
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